Resonant Two-Magnon Raman Scattering and Photoexcited States in 
Two-Dimensional Mott Insulators 
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We investigate the resonant two-magnon Raman scattering in two-dimensional (2D) Mott insula- 
tors by using a half- filled 2D Hubbard model in the strong coupling limit. By performing numerical 
diagonalization calculations for small clusters, we find that the Raman intensity is enhanced when 
the incoming photon energy is not near the optical absorption edge but well above it, being con- 
sistent with experimental data. The absence of resonance near the gap edge is associated with the 
presence of background spins, while photoexcited states for resonance are found to be character- 
ized by the charge degree of freedom. The resonance mechanism is different from those proposed 
previously. 

PACS numbers: 78.30.Hv, 71.10.Fd, 78.20.Bh 
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Two-dimensional (2D) insulating cuprates are typical 
Mott insulators. The understanding of their electronic 
states is of crucial importance since they are parent com- 
pounds of high-Tc cuprate superconductors. In the Mott 
insulators, low-energy physics associated with the spin 
degree of freedom is described by a 2D Heisenberg model. 
However, the nature of high-energy excited states across 
the Mott gap is less clear as compared with that of the 
low-energy spin states. 

Such excited states are usually investigated by opti- 
cal absorption measurements, where photons excite the 
system across the Mott gap. In addition to the optical 
absorption, resonant Raman scattering is useful for the 
investigation of the photoexcited states. In the scatter- 
ing, resonant effects are observed as enhancements of Ra- 
man scattering intensities associated with phonons and 
magnons when the incoming photon energy uj\ is tuned 
near the Mott gap. In the 2D insulating cuprates such as 
La2CuC>4, YBa2Cu306, and Sr 2 Cu02Cl2, the a>; depen- 
dence of two-magnon (2M) scattering intensity with B\ 
scattering geometry shows an interesting feature 0, |): 
there is no resonance at a peak near the gap edge in the 
absorption spectra (~1.8— 2 eV) but the resonance oc- 
curs at around io\ ~3 eV, well above the absorption-edge 
peak. Since such w\ dependence of the intensity is not 
a common behavior of transition metal oxides novel 
properties of the photoexcited states are expected in the 
2D Mott insulators. In contrast to a large number of 
theoretical studies about the asymmetric line shape of 
the B\ 2M Raman spectrum ||, [|, there are a smaller 
number of theoretical works on the cjj dependence of the 
intensity [t], |[ [|. As the origin of the resonance, a 
triple resonance mechanism based on a spin-density-wave 
formalism (?], ^] and an excitonic mechanism in which 
the charge-transfer (CT) exciton is taken into account ||] 
have been proposed so far. 

In this Letter, we propose a new mechanism of the 
resonant 2M Raman scattering and clarify the nature of 
the photoexcited states in 2D Mott insulators. Numer- 
ical calculations of the 2M Raman intensity for various 
scattering geometries as a function of u>i are performed 



for a half-filled 2D Hubbard model in the strong cou- 
pling limit. We find that the model can reproduce the 
experimental observations where the resonance in the B\ 
geometry occurs when uj\ is not near the Mott gap edge 
but well above it. The absence of resonance near the gap 
edge is found to be closely associated with the presence 
of background spins. On the other hand, the photoex- 
cited states where the resonance occurs show less impor- 
tance of the spin degree of freedom. Rather, the spatial 
distribution of charge carriers created by photons plays 
a dominant role in the resonance. It is crucial to no- 
tice that these conclusions are obtained by treating the 
Hubbard model exactly, taking into account full configu- 
rations associated with both the charge and spin degrees 
of freedom, in contrast to the previous works J?], |§, [§. 
Two-magnon Raman intensities with A\ , B2 , and A2 ge- 
ometries are found to show resonance behaviors that are 
different from those of B\. 

Insulating cuprates are known to be CT-type Mott 
insulators, where both Cu3<i and 02p orbitals partic- 
ipate in the electronic states. However, it is well es- 
tablished that the electronic states of the CT-type in- 
sulators can be described by a Hubbard model with a 
half- filled single band Jl(|. The Hubbard Hamiltonian 
with the nearest neighbor (NN) hopping in 2D is given 



by i?Hub = ~tJ2 



H.C.) + UY.i n iA n i,l, 



where c\ is the creation operator of an electron with 

spin cr at site i, ni t<r — c\ a Ci^, runs over pairs on 

the NN sites, t is the NN hopping integral, and U is the 
on-site Coulomb interaction. The value of t is estimated 
to be t ~ 0.35 eV (To). The value of U is estimated 
to be U =10t for the gap values to be consistent with 
experimental ones. 

In the strong coupling limit (U t), the ground state 
at half filling has one spin per site; i.e., there is no doubly 
occupied site. In this case, the low-energy excitation of 
the Hubbard model may be described by a Heisenberg 



model: H = J£ (lj) (S, ■ S 



iriiTij), where Si is the 
spin operator with S = 1/2 at site i, jii—ni^+ni.i, and 
J = At 2 /U. On the other hand, the photoexcited states 
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have both one doubly occupied site and one vacant site. 
In order to obtain an effective Hamiltonian describing 
the photoexcited states, we restrict the Hilbert spaces to 
a subspace with one doubly occupied site. Performing 
the second order perturbation with respect to the hop- 
ping term Ht in the Hubbard Hamiltonian, the effective 
Hamiltonian is given by 



H cS = TliHtlLx - -IIiZr t II 2 ff t IIi + -IlxHJIoHJLi + u, 

where Ho, III, and H2 are projection operators onto the 
Hilbert space with zero, one, and two doubly occupied 
sites, respectively. A complete expression of Eq. p) has 
been given elsewhere [ Q |. 

At zero temperature, the Raman scattering intensity 
is expressed in a Fermi's golden rule form, R(u>) = 
Ef \(f \ M n\ *)| 2 8 (w - Ef + Ei), where \ f) denotes a fi- 
nal state with energy Ef and |z) is an initial state with 
energy Ei. In the present case, \i) is just the ground 
state of the Heisenberg Hamiltonian Hq, while |/) is an 
excited state of Hq with two magnons. The Raman shift 
is given by u = cj; — u>{, u>{ being the outgoing photon 
energy. Mr is the Raman tensor operator. We first con- 
sider the Bi representation of the C± v group because its 
contribution dominates 2M Raman intensity in 2D Mott 
insulators ^ [l2|]. The dominant contribution to the B\ 
Raman spectrum comes from a matrix element, 



(/ 



Mi 



2 ^ 



(/ \jx I lA \jx \i)~{f \jy I H) (H \jy I i) 



En, - Ei - w; + iF 



(2) 



with the current operator j a = it • (c, 

cj a c i+a CT ). Here, the intermediate state is an eigen- 
state of H e g with the eigenenergy E^ and T is a damping 
rate for We note that the creation and annihilation 
operators, c\ a and c ia , in the current operator are pro- 
jected onto the subspaces with either zero or one doubly 
occupied site. 

We perform numerical calculations of the B\ Raman 
spectrum for a V20 x y20 square cluster with periodic 
boundary conditions. The Heisenberg ground state is cal- 
culated by the Lanczos method, and the spectrum R(to) is 
obtained by utilizing a modified version of the conjugate- 
gradient method together with the Lanczos technique. 

The inset of Fig. 1(a) shows a 2M Raman scattering 
spectrum with U/t=lO (J/i=0.4), where Wi=8i is chosen 
to be consistent with experimental conditions The 
spectrum is dominated by a large peak at w=1.37i=3.43 J 
accompanied by higher-energy small ones at to ~1.7i, 
being similar to that obtained by using a standard Ra- 
man spin operator 0] . The spectral weight changes with 
uii without notable change of line shape. The spectral 
weight of the u—1.37t peak, denoted as i?2mag, is plot- 
ted as a function of u>i in Fig. 1(a). i?2mag shows max- 
imum resonance at w;~llf. Comparing the w- x depen- 
dence with the optical absorption spectrum £2(^1), given 
by ll>^ 2 J2 u l(M Us I*) I ^(uJi-E^+Ei), we find that i? 2 ma g 
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FIG. 1: The dependence of (a) two-magnon Raman intensity 
J?2ma g with Bi geometry and (b) absorption spectrum £2 on 
incoming photon energy Ui in a 20-site cluster with U/t=W 
(J/t—OA). Inset of (a) shows two-magnon Raman spectrum 
at (vi=8t. i?2ma g represents the weight of the lowest-energy 
peak at u>—1.37t in the inset. The solid line in (b) is obtained 
by performing a Lorentzian broadening with a width of QAt 
on the delta functions denoted by vertical bars. 



is not enhanced when lui is tuned to the absorption-edge 
peak. Rather, the resonance occurs when oj\ is inside a 
continuum well above the Mott gap. The dependence 
of ^?2mag is thus similar to the experimental data men- 
tioned above Jl|, ^]. The energy difference between the 
absorption-edge and the resonance maximum is about 
4t~1.4 eV, being also consistent with experimental data 
with the difference of ~1.35 eV 0. We note that, in 
addition to the 20-site cluster, smaller clusters with 4x4 
and -\/T8x \/l8 sites show similar dependences (not shown 
here). We also checked that inclusion of longer-range 
hopping and Coulomb terms, which are used for a pre- 
cise description of cuprates, does not change the above 
results 
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The 2D Mott insulators before photoexcitation have 
localized spins interacting with each other via the antifer- 
romagnetic (AF) exchange interaction J. Charge carriers 
introduced into the insulators are known to induce a spin 
cloud around the carriers as a consequence of the mis- 
aligned spins along the carrier-hopping paths. It is thus 
natural to expect that the optical responses are strongly 
influenced by the presence of the localized spins. This 
can be checked by changing the value of the exchange 
interaction J/t. In Fig. 2, the absorption spectrum mul- 
tiplied by ujf (solid lines), a(u>i) = w ; 2 e2(u;i), is shown 
for various values of J, together with the 2M Raman in- 
tensity i? 2ma g (dotted lines). In the figure, the energy 
is shifted by U = At 2 / J. The absorption-edge peak is 
found to be very sensitive to the value of J: With de- 
creasing J, the peak intensity decreases, while the res- 
onant enhancement of i?2mag remains even for small J, 
keeping the energy difference between the resonance po- 
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FIG. 2: Absorption spectra a(u>i) and /3(u>i) from intital (solid 
lines) and final (dashed lines) states of the Raman process, re- 
spectively, together with two-magnon Raman intensity 7?2ma g 
(dotted lines) in a 20-site cluster for three J/t values. The 
incoming photon energy is shifted by corresponding value of 
U. The arrows denote positions of the absorption-edge peak 
and the maximum of ifemag- 



sition and the absorption edge almost constant. We note 
that such J dependence of the resonance is different from 
that expected from a triple resonance mechanism of the 
2M Raman scattering proposed based on a spin-density- 
wave (SDW) formalism for the half-filled 2D Hubbard 
model |?], H , where a maximum of resonance occurs near 
the upper edge of the absorption band with a width of 
8 J and thus the resonance energy must approach the 
absorption-edge with decreasing J. This implies that 
the origin of resonance we have found is different from 
the triple resonance seen in the SDW approach. 

The sensitivity of the absorption-edge peak to J seen in 
Fig. 2 is easily understood in the following way. The ini- 
tial state, i.e., the ground state of the Heisenberg model, 
is dominated by Neel-type spin configurations inducing 
AF order. The electric filed applied to the initial state 
creates a pair of empty and doubly occupied sites with 
one-lattice spacing, keeping the Neel-type spin configura- 
tions dominant in the spin background. The intensity of 
the absorption spectrum, which is given by | (/i \j x \i)\ , is 
thus predominantly governed by two factors in the pho- 
toexcited states (i) the probability of finding the 

two photocreated sites with one-lattice spacing, and (ii) 
the probability of the Neel-type configurations. In the 
photoexcited state at the absorption-edge peak |/x e d g e), 
the probability (i) is not largest among other possible 
two-site configurations with longer distances. Therefore, 
the Neel-type configurations in the spin background are 
crucial for the peak intensity. In fact, evaluating the 
wave function of |// e dge), we found that the probability 
(ii) increases with increasing J as expected. Therefore, 
at large J, the Neel configurations contribute to enhanc- 



ing the matrix element (^ dgo \jx \ i) , leading to the large 
absorption-edge peak. 

Using this consideration, it is easy to understand why 
the absorption-edge peak does not resonate with the 2M 
Raman scattering. In the Raman process, the matrix ele- 
ment of (/ \j x | fi) also contributes to the intensity, where 
|/) is the 2M final state. The spin correlation in |/) is dif- 
ferent from that in the initial state because two magnons 
are excited. By examining the wave function of |/), we 
found that the weight of the Neel-type spin configura- 
tions are very small in |/). Therefore, the matrix element 
(/ \jx\ A*edge) is very small, resulting in no resonance at 
the absorption edge. In order to confirm this, we plot in 
Fig. 2 an absorption spectrum from |/) (dashed lines), 

defined as /3( Wi ) = £ M |(/z \ Jx \ f)\ 2 Sfa - + Ek). As 
expected, it exhibits no enhancement at the photoex- 
cited state where a(u>i) shows the absorption-edge peak 
at fa - U)/t=-3A for J/*=0.5. 

In this context, recent resonant 2M Raman scatter- 
ing experiments performed on a ladder-type cuprate are 
interesting because the u>i dependence of the Raman in- 
tensity shows the same behavior as the absorption spec- 
trum |H| in contrast to the 2D systems. In the ladder 
system, the difference between the initial and final states 
is expected to be smaller than that in 2D because of a 
spin liquid ground state in the ladder. This may explain 
the difference between 2D and ladder fll6| . 

In contrast to the absorption-edge peak, a broad-peak 
structure, where the 2M Raman scattering resonates, is 
insensitive to J as shown in Fig. 2. This implies that 
the spin background has less effect on the structure and 
thus the relative position of the empty and doubly oc- 
cupied sites, i.e., the probability (ii) mentioned above, 
is important. In fact, by examining the two-site corre- 
lation with one-lattice spacing in all of the photoexcited 
states for a 4x4 cluster, we found that the correlation 
shows an enhancement around the broad peak structure 
(not shown here). The fact that 0{u>\) in Fig. 2 exhibits 
a similar broad-peak structure at the same energy re- 
gion, indicated by arrows, also implies little effect of the 
spin degree of freedom on the corresponding photoexcited 
states and thus supports the idea that the photoexcited 
states for the 2M Raman resonance are characterized by 
the charge degree of freedom. 

Finally, in order to clarify how the resonant behaviors 
mentioned above are influenced by scattering geometry, 
we examine 2M Raman spectrum for A\, B2, and A2 
symmetries. The corresponding Raman operators con- 
tain the current operators with the form of j x jx + jyjy, 
jxjy + jyjx, and j x j y - j y j x for the At, B 2 , and A 2 , re- 
spectively. The insets of Fig. 3 show the Raman spectra 
at LUi—8t. The intensity maximum in both the A\ and B2 
Raman spectra appears at w~1.7t, which is higher than 
that of B\ (1.37i) in Fig. 1. In addition, the spectra 
are broader than that of B\ . These features are qualita- 
tively similar to experimental data p"2[ . The A 2 Raman 
spectrum that reflects chiral spin fluctuations shows a 
very broad maximum at higher-energy u=3t, being also 
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FIG. 3: The dependence of integrated two-magnon Raman 
intensities up to u;=3.5t for (a) Ai, (b) B2, and (c) A2 geome- 
tries on incoming photon energy u>i in a 20-site cluster with 
U/t— 10 (J/t=0A). Insets show two-magnon Raman spec- 
trum at Ui—8t, where solid lines are obtained by performing 
a Lorentzian broadening with a width of OAt on the delta 
functions denoted by vertical bars. The dotted line in (a) 
represents integrated intensity of B\ up to u>=3.5t. 



consistent with the experiments ]12| . In the main panels, 
the oj\ dependence of the integrated Raman spectra up 
to u=3.5t is shown. For comparison, the integrated B\ 
intensity up to u>=3.5t is plotted in Fig. 3(a) as a dotted 
line. It is found that there is no pronounced resonance at 
the absorption edge of Wi=6.2i, being similar to the B\ 
geometry. On the other hand, the resonance behaviors 
at higher lo\ are different from that of B\ : the resonance 
maxima do not appear at Wi~10i and the intensities are 
smaller than that of B\ geometry. Systematic experi- 
mental studies of resonant 2M Raman scattering decom- 
posed into various scattering geometries for a wide range 
of LJi are necessary in order to confirm these differences. 



We note that, from detailed analyses that are similar to 
those done for the B\ geometry, the main origin for the 
resonance maxima is found to be the same as that of 
Bi, i.e., the charge degree of freedom. For example, a 
resonance maximum at tJi~12i in A\ remains even for 
small values of J and the energy of the maximum ap- 
proaches the center of the photoexcited band = U), 
where the probability of finding the two photocreated 
sites with one-lattice spacing is very large. 

In summary, we have clarified the mechanism of the 
resonant two-magnon Raman scattering and the nature 
of the photoexcited states in 2D Mott insulators. Numer- 
ical calculations for a half-filled 2D Hubbard model in the 
strong coupling limit have clearly reproduced the exper- 
imental observations that the resonance occurs when the 
incoming photon energy is well above the Mott gap en- 
ergy, but does not when the energy is at the absorption- 
edge peak. The spin degree of freedom is found to play 
a crucial role in the absence of resonance at the peak. In 
particular, the difference of spin configurations between 
the initial and final states is a main reason for the absence 
resonance. In contrast, the photoexcited states where the 
resonance occurs show less influence of the spin degree of 
freedom. Rather, the spatial distribution of charge car- 
riers created by photons plays a dominant role in the 
resonance. In other words, the charge degree of free- 
dom controls the resonance of the two-magnon Raman 
scattering in the 2D Mott insulators. The mechanism of 
the resonance and nature of the photoexcited states that 
we have found are different from those proposed previ- 
ously @ Hi ||]- It should be emphasized that, in contrast 
with the previous works, our conclusions are obtained by 
treating the Hubbard model exactly, taking into account 
full configurations associated with both the charge and 
spin degrees of freedom. 
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